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Abstract 

We investigate scalar and tensor fields in the brane model solution for the 5D space- 
time with standing gravitational waves in the bulk and show that even in the case of 
increasing warp factor there exist normalizable zero modes localized on the brane. 

PACS numbers: 04.50.-h, ll.25.-w, 11.27.-Fd 
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1 Introduction 



The scenario where our world is associated with a brane embedded in a higher dimensional 
space-time with non-factorizable geometry [T],|2] has attracted a lot of interest with the aim of 
solving several open questions in modern physics. Most of these models were realized as time 
independent field configurations. However, there have appeared several braneworld models that 
assumed time-dependent metrics and fields |3J. Here we consider the braneworld scenario with 
non-stationary metric coefficients recently proposed in [4]. The braneworld is generated by 5D 
standing gravitational waves coupled to a phantom-like scalar field in the bulk. 

A key requirement for realizing the braneworld idea is that the various bulk fields be localized 
on the brane. For reasons of economy and avoidance of charge universality obstruction [S] 
one would like to have a universal gravitational trapping mechanism for all fields. However, 
there are difficulties to realize such mechanism with exponentially warped space-times. In the 
existing (l+4)-dimensional models spin and spin 2 fields can be localized on the brane with 
the decreasing warp factor |2], spin 1/2 field can be localized with the increasing factor [6], and 
spin 1 fields are not localized at all [7]. For the case of (l+5)-dimensions it was found that spin 
0, spin 1 and spin 2 fields are localized on the brane with the decreasing warp factor and spin 
1/2 fields again are localized with the increasing factor [S]. There exist also 6D models with 
non-exponential warp factors that provide gravitational localization of all kind of bulk fields on 
the brane p|, however, these models require introduction of unnatural sources. 

As it was mentioned in [1] standing wave braneworld model can provide universal gravi- 
tational trapping of zero modes of all kind of fields in case of rapid oscillations of standing 
gravitational waves in the bulk. To understand the idea we want to remind that standing 
electromagnetic waves, so-called optical lattices, can provide trapping of various particles by 
scattering and dipole forces [10]. In [11] localization was demonstrated through quadruple 
forces as well. It is known that the motion of test particles in the field of a gravitational wave 
is similar to the motion of charged particles in the field of an electromagnetic wave [12] . Thus 
standing gravitational waves could also provide confinement of matter via quadruple forces. As 
a classical example let us consider the equations of motion of the system of spinless particles 
in the quadruple approximation [13] . 



where is the total momentum of the matter field and J"'/^^^ is the quadruple moment of 
the stress-energy tensor for the matter field. The oscillating metric due to gravitational waves 
should induce a quadruple moment in the matter fields. If the induced quadruple moment is 
out of phase with the gravitational wave the system energy increases in comparison with the 
resonant case and the fields/particles will feel a quadruple force, -F^^^, which ejects them out 
of the high curvature region, i.e. it would localize them at the nodes. 

In this Latter we show that standing gravitational waves induced by the brane in the 
bulk really can localize scalar and tensor field zero modes on the brane even for the case of 
exponentially increasing warp factor. 



Dp'' _ 
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2 The model 



In this section we briefly recall the 5D standing wave braneworld model proposed in The 
braneworld is generated by gravity coupled to a non-self interacting scalar phantom-like field 
which depends on time and propagates in the bulk. The action of the model has the form: 



(2) 



where and A5 are 5D Newton and cosmological constants respectively. 

To avoid the well-known problems with stability, which occur with ghost fields, we can 
associate the ghost-like field with the geometrical scalar field in a 5D integrable Weyl model 
|15j . In the Weyl model a massless scalar appears through the definition of the covariant 
derivative of the metric tensor, 

DaQmn = gMNdA4> ■ (3) 

This is a generalization of the Riemannian case where the covariant derivative of the metric is 
zero. The assumption implies that the length of a vector is altered by parallel transport. 
Weyl's scalar field in ([3]) may imitate a massless scalar field - either an ordinary scalar or 
ghost-like scalar. The gravitational action for the Weyl 5D integrable model can be written as 



S = J d^'x^ 



^ {R - 2A5) - (6 - ^i)g''''dM(pdN^ 



(4) 



where ^ is an arbitrary constant. For example, for = 11/10 the action (jlj) exactly coincides 
with ([2]). Thus we can start with a 5D Weyl model and require that we have Riemann geometry 
on the brane by assuming that the geometrical scalar field (j) is independent of 4D spatial 
coordinates, x*, and vanishes on the brane. The definition of via ([3]) avoids the usual instability 
problems of ghost fields since the geometrical fields have specific couplings with matter fields 
and it is known that the Weyl model is stable for any value of ^. 
The Einstein equations for the action (|2]) read: 

Rmn — -^QmnR = ^ttG^Tmn — ^5gMN ■ (5) 

The metric ansatz we use is: 

ds^ = e^'^l"! {dt^ - e"c/x2 - e'^dy^ - e'^'^dz'^) - dr^ , (6) 

where a is a constant, which corresponds to brane width. The peculiarity of the model (jS]) is 
that the brane, located at r = 0, possesses anisotropic oscillations and send a wave into the 
bulk (as in [16]), i.e. the brane is warped along the spatial coordinates x, y, z through the factor 
^u{t,r) ^ which depends on time t and the extra coordinate r. 

On the other hand, the phantom-like scalar field r) obeys the Klein-Gordon equation 
on the background space-time given by (EI), 

^ 5m(v^/'^9^0) = e-^'^l^l^ - 0" - 4a6(r)0' = , (7) 
V 9 
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where e(r) is the sign function and overdots and primes mean derivatives with respect to t and 
r respectively. 

We further rewrite the Einstein equations in the form: 

2 

Rmn = -OmuOno- + -gMN^5 , (8) 
where the gravitational constant has been absorbed in the definition of the scalar field: 

a = v/SttGs . (9) 
It turns out that the system is consistent only when the fields a and u are related [1]: 



u{t, r) 



(10) 



The anisotropy of the metric ([6]) forces us to define the energy-momentum tensor of the 
brane with different tensions along different directions: 



r^' = S{r) diag[Ai, A^, Ay, A^, 0] , (A^ = \y 



'ID 



The tensions Am (m = t, x, y, z) in general depend only on time. Taking into account (|TT]) the 
field equations ([8]) change as follows: 



Rmn = -dM(ydN(y + -^gMN^b + 8'kG5Tmn 



where the reduced energy-momentum tensor, 
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tmn — tmn — -^Qmn't , 
corresponds to the matter content on the brane and takes the form: 

TMN = ^Sir) diag [2\t + 2e-"A, + e^^A,, e^A* + A, - e^^A, 
e"At + A, - e=^"A„ e'^^At - 2e-=^"A, + 2A„ O] 
The non-zero components of the Ricci tensor for the metric (EI) read: 



Rtt 




Rxx 


= Ryy = 


Rzz 


= g2a|r|- 




= --u" 
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3 . 



^2a\r\+u 



^e-^'^l'^lu - - 2ae{r)u' - 2a6{r) - ^u" 



ir)+u"] 



(12) 



(13) 



(14) 



(15) 



Using the fine tuning: 

As = , (16) 
the system of Einstein equations (1121) reduces to a single ordinary differential equation, 

u-u" - Aae{r)u' = . (17) 
In addition, for the brane tensions we have [17J: 

u{t,0) / 1 \ 

A. = ^^ = ^(6a--K]j5(r) , (18) 

-2u{t,0) 

K = ^-7^(6a+K])5(r) , 
A, = 0, 

where u (t, 0) and [u'] denote the value of u{t,r) and the jump of its first derivative at r = 0, 
respectively. 

The standing wave solution to f|T7|) can be constructed by implementing the ansatz : 

M(t,r) = Csin(a;t)/(r) , (19) 
where C and u are real constants. From f|T7|) we get the equation for the radial function: 

/" + 4ae(r)/' + u^'^''^''^ f = . (20) 
The general solution to this equation has the following form: 



/(r) = e 



-2a\r\ 



A J2 I — e-"^!^! 1 + 5 F2 
a 



(21) 



where A, B are arbitrary constants and J2 and Y2 are second-order Bessel functions of the first 
and second kind, respectively. 

As pointed out in along with the solution fl2T]) . the ghost-like field cr(t, r) must be 
unobservable on the brane. Taking into account the relations ( ITOl) and (1191) we can accomplish 
this requirement by setting the boundary condition: 

/WUo = 0. (22) 

Since J2 and Y2 are oscillatory functions, in the case of increasing (decreasing) warp factor for 
some fixed values of the constants A and B the function /(r) can have finite (infinite) number 
of zeros. Thus the above boundary condition can be written in the form which quantizes the 
oscillation frequency, of the standing wave in terms of the brane width, 1/a, i.e. 

V-,=Xn, (23) 
\a\ 

where X„ is the n^^ zero of the function /(r). Correspondingly, the functions cr(t, r) and u{t^ r) 
vanish at the finite (infinite) number of points along the extra dimension. These points - 
the nodes of standing wave - can be considered as 4D space-time 'islands', where the matter 
particles are assumed to be bound. 
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3 Localization of scalar and metric fields 



Since the equations for extra components of tensor and scalar fluctuations are similar for the 
simplicity let us consider localization problem only for scalar fields, which are defined by the 
5D action: 

S^ = -\j Vgdx'dr g^'^dM<l>dN<^ . (24) 
Corresponding Klein- Gordon equation has the form: 

^ Om (v/^/^^9^$) = . (25) 

Determinant for our ansatz ([6]) is equal to = e^"'^' and fl25|) can be written as: 

[a^ - e-"(a2 + dl) - e^"^^] $ = e^'^l^l (e^"l''l<l>')' , (26) 
where primes denote derivatives with respect to r. Let us look for the solution in the form: 

$ (t, X, y, z, r) = m{t, r)x{x, y)i{z) , (27) 
which transforms fl26|) into the system of the equations, 

{dl + d'y)x+{pl+pl)x = 0, 

dl^ + p'A = 0, (28) 
d^-^ + [{pI + p2)e-" + py^] ^ = e^'^l'-l (e^'^l^'l^')' . 

In the case u ~ the parameters Px,Py and Pz can be regarded as momentum components along 
the brane. Their exact physical meaning is not clear in general case since fl28|) are equations 
with variable coefficients. 

To solve the system ( l28l) we assume that the constant A in the solution ( 12T]) is zero and the 
oscillation function of gravitational waves has the form: 

u {t, r) = B sin {ut) e^^^I'lFa e"'^!''!^ , (29) 

where i? is a constant and Y2 is the second-order Bessel function of the second kind. Physical 
solutions described by Bessel functions usually contain the first kind functions, J„, because 
of their regularity at the origin. However, in our case we can use even Bessel function of the 
second kind, since the argument in f l29|) . e~"l^la;/|a|, is always positive and becomes zero 
only for the case of increasing warp factor (a > 0) at r — t- ±00. 

To show explicitly localization by standing waves we consider the extreme case with the 
increasing warp factor a > 0. The function (129|) is zero at the position of the brane, r = 0, due 
to fine tuning ( |23|) . where X„ now is one of the zeros of 1^2- 

For 4D massless scalar modes on the brane we have the following dispersion relation: 

E''=pl+pl+pi. (30) 
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When the frequency of standing gravitational waves is much larger than frequencies associated 
with the energies of the particles on the brane, 



(31) 



we can perform time averaging of oscillating exponents in the equation for scalar field fl28|) . 
Using the expressions: 



2-k/uj 



2^ 







[sin(ci;t)]™ dt 



2-2" (2n)! (n!) 



(m = 2n + 1) 
(m = 2n) 



bu 



E 

n=0 



{buy 



(32) 



where h is an integer, for the time averages of oscillating exponents f l32|) we get: 



(bBf^e 



2n„—4an\r\ 



n=l 



( ^e- 

a 



-a\r\ 



2n 



(33) 



where Jq is the modified Bessel function of the zero order with the argument: 



z = \bB\ 

Now separating variables, 



2a 



— e 



(34) 



(35) 
(36) 



for the extra factor p(r) of the wave function we obtain: 

(e4«Mp')'_e2«MF(r)p = , 

where 

F(r)=((e->-l)(p^+p^ + ((e->-l)p^ 

It is more convenient to put Eq. (!35|) into the form of an analogue non-relativistic quantum 
mechanical problem by making the change: 

p(r) = e-2'^l'''^/'(r) . (37) 

For ilj{r) we find: 

^" - U{r)^lj = , (38) 

where the function, 

f/(r) = 4a(5(r) + Aa^ + e-^^^'^ F (r) , (39) 

is the analog non-relativistic potential. FIG. 1 shows the behaviour of U (r) in the case of first 
zero of the Bessel function Y2, i.e. when 



- ^ 3.38 . 

a 



(40) 



12 n 



2- 



I 1 1 ' 1 ' 1 ' 1 ' 1 ' 1 

-3-2-10 1 2 3 

r 

Figure 1: The total gravitational potential for scalar particles in the bulk, fl39l) . 



To study general behaviour of extra part of the scalar zero mode wave function we explore 
Eq. f l38|) in two limiting regions, far from and close to the brane. 

Far from the brane, r — )■ ±oo, Eq. f l38|) has the following asymptotic form: 

^" - 40^^ = , (41) 

with the solutions if) ~ e^^^'^L Taking into account the normalization condition for the wave 
function and a > 0, we choose if) ~ e~^"'^'. So, far from the brane for the extra factor of the 
wave function ( 1371) we obtain: 

p(r).-.±oo ~ e-^'^l'-l . (42) 
In the second limiting region, r — )■ ±0, i.e. near the brane, Eq. fl38|) can be approximated 

as: 

ip" - [Aa5{r) + 40^ + cV^] ^ = Q , (43) 
where we have introduced the constant: 



The solution of 



is: 



if) (r) 



D 



W , -, c r 



c '4' 



Yd- 

a 



+ 



1 2 



(44) 



20F 



2v^ 



r(f + f) V-c^ 



M(--,-,c|r| 



(45) 
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where Z) is a constant, T is the gamma function and W and M denote Whittaker functions. 
Imposing boundary condition at the origin, 



P\r)\ 



r=0 







close to the brane extra part of scalar wave function will behave as (see FIG. 2): 



r(l + ^) 



1 + 



T2~ 



91 I ^ 

ac \r\ 



15 



+ O (r^) 



(46) 



(47) 



So the extra part of scalar field zero-mode wave function p(r) has maximum at r = 0, falls 
off from the brane and turns into the asymptotic form (H2|) at the infinity. 
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Figure 2: The extra part of scalar zero mode wave function close to the brane, (H71) . 



In the action of scalar fields flM|) the determinant, ^/g = e^"^'^', and the metric tensor with 
upper indices, g^^^ , give the total exponential factor e^'^''"', which obviously increases for a > 0. 
This is the reason why in the original brane models |2] the scalar field zero modes, with the 
constant r-depended extra part of the wave function, can be localized on the brane only in the 
case of decreasing warp factor (i.e. a < 0). In our model the extra part of wave function fl57|) 
is not a constant, moreover, for a > it contains the exponentially decreasing factor e"^"'*^'. 



4 Conclusion 

In this Latter we have investigated localization problem of scalar and tensor fields within the 5D 
standing wave braneworld proposed in [1]. Since equations for the extra part of wave functions 
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of scalar and tensor fields are similar for simplicity we have explicitly considered in details only 
the case of scalar fields. We have shown existence of normalizable zero modes on the brane 
even for the case of increasing warp factor. Further investigations towards the localization of 
spinor and vector fields in this model are in progress. 

Acknowledgments: The research was supported by the grant of Shota Rustaveli National 
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